Modified dynamics of weakly coupled BEC's Josephson Junction (BJJ) 
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The tunnelling quantum dynamics of bimodal BJJ system is modified through introducing an 
equilibrium condition, which is based on the assumption that the BJJ is tend to keep on its ground 
state (with a lowest energy) during the oscillation. The tunnelling dynamics of BJJ with symmetric 
and asymmetric traps is discussed through numerically solving the modified equations. Station- 
ary states are found to exist in the both BJJs. Compared to previous works, the macroscopic 
quantum self trapping (MQST) is auto-avoided. Meanwhile, it is revealed that the BJJ oscillates 
with its inherent frequency which is only related to the Josephson energy, which has been testified 
experimentally in other contexts. 
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o 



i 

C 

o 
o 



> 

OS 

o 
o 

On 
O 

m 
o 



T3 
C 

o 
o 



I. INTRODUCTION 

Since Bose-Einstein condensation was first detected in 
1995 0, there have been rapid and critical developments 
in experimental techniques 0, 0, 0, IE ■ In 1997, in- 
terference fringes in two overlapping condensates were 
observed^. Consequently in 1998, the superposition of 
condensed atoms in different hyperfine levels was cre- 
ated 0, And the evolution of the relative phase of 
two coupled condensates was measured by interferome- 
try techniques 0,0- I n 2001, the direct observation of 
an oscillating atomic current in one-dimensional array 
Josephson junctions was realized With these precise 
manipulation of Bose-Einstein condensates (BEC's), it 
has enhanced the possibility of tailoring the new quan- 
tum systems to a degree not possible with other quan- 
tum systems, like superfluid. And the studies of spatial 
coherence naturally raise the question of measurement 
and exploration of temporal phase coherence between two 
condensates. 

Theoretically, with the approximation of non- 
interacting atoms and small-amplitude Josephson oscil- 
lations, some aspects of temporal phase coherence have 
been investigated in the context of BEC by means of 
Josephson junctions 0, 0- Yet many important fea- 
tures of this subject remained to be explored. A non- 
excitation system of interacting bosons confined within 
an external potential can be described by a macroscopic 
wave function, which has the meaning of an order pa- 
rameter and satisfies the nonlinear Gross-Pitaevskii (GP) 
equation [Toj . The semiclassical tunnelling quantum dy- 
namics using GP equation has been studied |ll|. And 
many interesting phenomena such as macroscopic quan- 
tum self-trapping (MQST) [H [H and tt states 
[TTL fl3| were predicted in a weakly coupled double BEC. 
In those works, the discrete nonlinear Schrodinger equa- 
tion (DNLSE) has served as a powerful tool for describ- 
ing the boson Josephson junction (BJJ) jUEJ. But so 
far, many investigations have cast doubts on the validity 
of the DNLSE for conserved quasiparticles interacting 
with a boson field [l4|- The objections are not to the 



GP equation itself , but to the quasiparticle approxima- 
tion instead. The earlier results at the condition of short 
time scales obtained through the DNLSE in the context 
of coupled quasiparticle-boson systems were verified by 
the fully quantum version of the BJJ tunnelling model 
0,13- But on long time scales, the MQST in symmetric 
traps, which is predicted by DNLSE, was found to be 
destroyed in fully quantum dynamics |fj. 

As revealed in other context 0, the MQST (or z- 
symmetry breaking) in a translationally invariant Hamil- 
tonian is an artifact of the factorization assumption in- 
herent in the semiclassical dynamics. In the DNLSE pro- 
cedure, atoms are assumed to be localized in either of 
the traps with discrete wavefunctions. Then, the quan- 
tum states for the total BJJ system is the superposition 
of those eigenstates in each trap individually. Further- 
more, the population distribution is determined through 
the calculation of time dependent Schrodinger equation. 
Compared to the fully quantum dynamics, in the semi- 
classical description, the quantum thermodynamic equi- 
librium for BJJ system is overlooked. And a stationary 
BJJ is not restricted to its ground state with a lowest en- 
ergy. A consequent result is the artifact of MQST, which 
would not break down in the semiclassical description 
even on the condition of long time scales. 

Compared to previous works, the BJJ is assumed to 
occupy its ground state in present study. And a cor- 
responding lowest energy condition was deduced, based 
on which the semiclassical dynamics was modified. The 
numerical calculation of the modified dynamics revealed 
that the MQST was auto-avoided. Also, it was found 
that the BJJ oscillated with an inherent frequency, which 
was only related to the Josephson energy. 

In the following Sec. ^ the equilibrium behavior of 
BJJ is discussed, and the lowest energy state is found 
under a condition of a fixed relationship between the pop- 
ulation distribution and the quantum states of BEC's. In 
Sec. IIIII the semiclassical dynamics is modified by intro- 
ducing the stationary condition into the DNLSE proce- 
dure. In Sec. IIVI the tunnelling dynamics for BJJ with 
symmetric and asymmetric traps is discussed and some 
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numerical results are given. 



II. THE EQUILIBRIUM BEHAVIOR OF 
WEAKLY COUPLED BJJ 



The effective many-body Hamiltonian describing 
atomic BEC in a double-well trapping potential Vt rn.n{ r) 
can be written in the second-quantization form as 



H = J d 3 r*t[_|lv 2 + y irap ]* + ^ J d 3 r*t*t>H,.(l) 

Here m is the atomic mass and Uq = 4:nah 2 /m (a denotes 
the s-wave scattering length, measuring the strength 
of the two-body atomic interaction). The Heisenberg 
atomic field operators W and ^> satisfy the standard 
bosonic commutation relation [^'(r,t),^{r',t)] — S(r — 
r'). In the quasiparticle approximation, one can expand 
the field operators in terms of two local modes 



(2) 



Here fit (i)] = and ipi(r) stands for the local 



mode function of either well and satisfies 



d A ri)*{r)ij)i(r) w 



(3) 



Substituting Eq. @ into Hamiltonian QJ, the two- mode 
approximation of H is yielded 



H= ( E i & l^ + Kajaja^) - (Ja\a 2 + J*fiifi|), ( 4 ) 

i=l,2 

where the parameters are estimated by |l2j 

= [ efV^hl-V 2 + V trap ]A, 
J 2m 



U 



d 3 r\iPi\ 



(5) 



<Pr[ — Vlpl ■ Vlp 2 + VtraplptH 
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Here the interactions between atoms in different wells arc 
neglected as in the weakly coupled BEC's. Defining two 
local number operators hi — fitfi*, it is easy to verify that 
the total number operator N = fijfii + fijfia — Nt rep- 
resents a conserved quantity. After neglecting the con- 
stant term, the two-mode Hamiltonian H can be rewrit- 
ten as Ha: 



H' = E c (h - h g f - ( Jfi*fi 2 + .Phial), 
h = {fix - h 2 )/2, 

E c = \ 1 +\ 2 , (6) 

% = - ^i°) + (Nt - l)(Ar - A 2 )]. 

Here h is the number difference operator, E c is the 
"charging energy" and h g is known as the "gate charge" 

M 

Generally, the BJJ system woks in Fock regime, which 
gives E c 3> J, and Hamiltonian can be well approxi- 
mated by 



H = E c {h — hg 



(7) 



A BJJ system in ground state should yield a minimized 
energy (or say a maximized entropy) and the Hamilto- 
nian should be approximate zero, which gives 

h = h g . (8) 

Eq. (JSJ denotes the relationship between the population 
distribution and "gate charge" of an stationary BJJ sys- 
tem in its energy ground state. It is important to notice 
that the "gate charge" , which is related to the eigenstates 
of BEC in the either trap, is determined by two compo- 
nents, i.e. the boson-field interaction and self interaction, 
as shown in Eq. 10. In the present work, the number 
difference h = h g is not necessarily an integer, and the 
artifact of the factorization of assumption inherent in the 
semiclassical dynamics is automatically avoided. In this 
work, the boson-field interaction energy is considered as 
independent of the number difference n, while the self in- 
teraction energy is shifted by the atom numbers in cither 
trap [Hj . 



III. MODIFIED TUNNELLING DYNAMICS 

In the DNLSE procedure, the wavefunction Eq. © is 
rewritten in the quasiparticle form, with 



i=l,2 



(9) 



where N and 9 denote the population and phase shift 
of BJJ, and 4>i represents the local mode function fij(i). 
Then, the Hamiltonian © is reformed to 



i fi- 



at 



(E° + XiNi)<fn - Jfo 



ih-^ = {El + \ 2 Ni)4> 2 - J<f>i. 



(10) 
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Defining the fractional population imbalance and rela- 
tive phase as 

z(t) = (Nxit) - N 2 (t))/N T , 

/3(i)=0i(i)-0 2 (i). (11) 

Eq.(lO) becomes 

z(t) = -v/l-z(i) 2 sin[/?(i)], 
Kt) = AE/(2J) + jM= cos[/?(i)]. (12) 

Here t has been rescaled to a dimensionless time t2J/h. 
And A.E denotes the particle energy difference between 
atoms in either traps, which acts as an asymmetric pa- 
rameter in the tunnelling quantum dynamic of BJJ. In 
present work, the BJJ system is assumed to stay at its 
ground state. So the stationary condition Eq. © should 
be met. By introducing Eq. JSJl, Eq. (|T^|l should be 
modified. It was found that the dynamic equations keep 
its form as shown by Eq. (|12|1 . but with the energy dif- 
ference modified to 



where z(0) and /3(0) are the initial values of z and (3 
respectively. In the earlier work ^l| , the energy gap AE 
is the function of population imbalance z , which would 
result in some certain initial value of z satisfying Eq. 
(|16fl . Thus the BJJ system would maintain in MQST 
state even at the condition of long time scales. Compared 
to the previous work, a equilibrium condition Eq. © was 
introduced to modify DNLSE procedure in present work. 
It was found the artifact induced by DNLSE is auto- 
avoided because in the BJJ with symmetric traps the 
expression Hq = — \J\ — z cos (3 < 1 is satisfied by any 
given z(0) and /3(0). Naturally, because the tunnelling 
quantum dynamics behavior is restricted by equilibrium 
condition, it is a self-consisted result of Eq. JBJ) and Eq. 
()12J) that the oscillation frequency should be independent 
of initial condition z(0) and /3(0). 

IV. TUNNELING DYNAMICS OF BJJ WITH 
SYMMETRIC AND ASYMMETRIC TRAPS 

A. The symmetric trap case 

For a symmetric BEC Josephson Junction (BJJ), the 
motion equations Ijl2(l is given by 



AE = (El - E$)/2. (13) 
And the total conserved energy is given by 



H = -^-AE-z- z 2 cos f3. (14) 

Eq. (|12|l can be rewritten in the Hamiltonian form ^l| 



dz 
dt 



dt 



dH 

dH 

dz 



(15) 



In deducing Eq. I|15l) . we have put forward the fact 
that at low oscillation frequency, the BJJ system would 
stay in its lowest energy state, say ground state, and the 
population distribution would be related to the "gate 
charge" . In the previous tunnelling quantum dynamic 
description of BJJ, the GP equation describing the mean- 
field dynamics of a BEC is reformed to the bimodal dis- 
crete nonlinear Schrodinger equation, and the calculation 
is carried out by solving self-consistent equation in each 
trap respectively 01 ■ ^ would bring about artifact since 
the equilibrium behavior has been overlooked. And a 
consequent result of previous tunnelling dynamics is the 
MQST 0. Generally, MQST is based on the condition 
that the initial conserved energy 



m = 



yi-z(t) 2 sin[/3(t)], 
z(t) 



:COB[/3(t)], 



(17) 



with a conserved energy H = — \/l — z 2 cos (3. The 
ground state is a symmetric stationary solution of Eq. 
(|17fl . with a conserved energy and i?_ = — 1 and 



2n7r, 
z = 0. 



The next stationary state with higher energy E. 
is an antisymmetric solution with 



[2n + 1)tt, 
z = 0. 



(18) 
= 1 



(19) 



The imbalanced initial population, i.e. z(0) ^ 0, would 
result in an oscillating solution of Eq. I|17|) , as shown in 
Fig. ^ It is seen that for a given z(Q), the popula- 
tion imbalance oscillates around it's zero-value, and the 
z-symmetry of bi-mode BJJ is preserved. Meanwhile, the 
oscillation amplitude A is determined by the initial pop- 
ulation imbalance z(0), i.e. A — \z(0)\. 



B. The asymmetric trap case 



H o = H[z(p),0(p)] > 1, 



(16) 



When the trap is asymmetric, the energy difference is 
not equal to zero, i.e. AE ^ 0. The numerical solution of 
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FIG. 1: Population imbalance as a function of dimensionless 
time t2J/h. The BJJ is of symmetric traps with initial relative 
phase /3(0) = 0, and the initial population imbalance z(0) 
takes the values 0, 0.1, 0.4 and 0.7 respectively. 

Eq. l|12fl is given in Fig. 2 with the trap asymmetry num- 
ber AE/2J = 1 and z(0) = 0, 0.4, and 0.7 respectively. 
It is shown that the population imbalance z(t) oscillates 
around its mean value z with an amplitude A — \z(0) — z\. 
Also the z(t) plotted against relative phase (3(t) is given 
in Fig. 3. It is expected that when z(0) = z = 0.447, 
oscillating amplitude would be zero, i.e. A = 0, and the 
BJJ system should arrive at a stationary state. The cor- 
responding ground state is obtained with the symmetric 
phases, 

p = 2mr, 

z = z. (20) 
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FIG. 2: Population imbalance as a function of dimensionless 
time t2J/h. The BJJ is of asymmetric traps with AJ5/2J = 
—0.5, and a initial relative phase /3(0) = 0, with the initial 
population imbalance z(0) takes the value 0.1, 0.4, 0.447and 
0.7 respectively. 



FIG. 3: Population imbalance as a function of relative phase 
Pit). The BJJ is of asymmetric traps with AE/2J — —0.5, 
and a initial relative phase /3(0) = 0. As the initial population 
imbalance z(0) approach a mean value of about z — 0.45, the 
oscillating amplitude approach zero. 

The next stationary state with higher conserved energy 
is given by 

(3= (2n+l)7r, 

z = z, (21) 

with antisymmetric phases. Moreover, as shown in Eq. 
(|12fl . the stationary state may be arrived only on the 
condition of z(0) = I, and is independent of its initial 
relative phases. 

It is important to notice that, both in Fig. ^and Fig. 
12 the oscillation frequency / is independent of either the 
trap asymmetry number AE/2J or the initial value z(0). 
It is the inherent character of BJJ, which is only related 
to the Josephson energy J, with 

/ cx - (22) 

Since Josephson energy is determined by the interwell 
potential height, it is experimentally feasible to testify 
the modified dynamics through Eq[221 In fact, it has pro- 
vided a testification of the present scheme [l^. 

V. CONCLUSIONS 

In conclusion, the tunnelling quantum dynamics of bi- 
modal BJJ system is modified under an assumption that 
the BEC's is tend to keep on the ground state with a 
lowest energy. And the lowest energy is determined by 
the population distribution. The dynamics of BJJ with 
symmetric and asymmetric traps is discussed with the 
modified equations. Compared with previous works, the 
MQST is auto-avoided in this work, as has been testified 
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by previous experimental works. Also, it is revealed that 
the BJJ oscillates with its inherent frequency, which is 
only determined by the Josephson energy. This has been 
testified experimentally. To the end, it is necessary to 
point out that it remains to be seen how quasipaticle ex- 
citation and energy dissipation of BEC would affect the 
tunnelling dynamics. 
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